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ABSTRACT 


The  on-line  operation  of  a  Fisher  detector  has  been  evaTrsted 
at  the  Tonto  Forest  Observatory  short-period  seisi  v  array.  :s  ng 
31  seismometers  and  one  fixed  threshold,  the  90  percent  incremental 
detection  threshold  on  the  Kuril  Islands  beam,  centered  at  A  =  70°, 
is  m^  =  4.15  +_  .05;  with  an  experimentally  determined  false  alarm 
rate  of  0.17  per  day.  This  compares  favorably  with  the  capabilities 
of  a  human  operator.  Storms  in  the  Kurils  significantly  affect  the 
distribution  of  amplitudes  of  the  Fisher  detection  trace;  and  we 
estimate  that  most  of  the  false  alarms  observed  at  the  operating 
threshold  can  be  traced  to  the  statistical  bias  introduced  by 
this  storm-generated  energy.  If  the  threshold  were  ad'usted 
to  maintain  a  constant  false  alarm  rate,  the  maximum  effect  on  the 
threshold  magnitude  would  be  0.05  m,  . 


TABLE  OF  CONTENTS 


Page  No 

ABSTRACT 

INTRODUCTION  1 

THEORY  6 

PROCEDURES  n 

RESULTS  14 

DISCUSSION  19 

ACKNOWLEDGEMENTS  21 

REFERENCES  22 

APPENDICES 
APPENDIX  I 

Calculation  of  the  F- Dis tribution 
APPENDIX  II 

Operating  Schedules  and  Parameters 


I  1ST  01;  FIGURES 


Figure  Title 

F  vaiues  as  a  function  of  false  alarm  rate  in 
false  alarms  per  day  lor  arrays  with  different 
numbers  of  elements. 


Figure  No 


False  aiarms  per  day  as  a  function  of  probability 
of  detection  for  certain  fixed  signal  to  noise  ratios 

U  the  nrray  3I?d  f°r  3  de^rees  of  freedom 

in  the  numerator,  3  (N-lj  in  the  denominator  where 

N  is  the  number  of  elements  in  the  array.  The' 

igure  is  for  an  array  of  3  elements. 

Raise  aiarms  per  day  as  a  function  of  probability 
of  detection  for  certain  fixed  signal  to  noise  ratios 
the  array  beam,  and  for  3  degrees  of  freedom 

i  ?h:UIlat0r’r3(f'1]  in  the  denominator  where 
is  the  number  of  elements  in  the  array.  The  figure 
is  for  an  array  of  5  elements. 

of1dete^I’mS  r  Jay  35  3  function  of  probability 
of  detection  for  certain  fixed  signal  to  noise  ratios 

on  the  array  beam,  and  for  3  degrees  of  freedom 
in  the  numerator  3  (N-l)  in  the  denominator  where 
is  the  number  of  elements  in  the  array.  The  fieure 
is  for  an  array  of  7  elements. 

of!6®!;!™5  ?er  da>’  as  a  function  of  probability 
o  detection  for  certain  fixed  signal  to  noise  ratine: 
on  the  array  beam,  and  for  3  degrees  of  freedom 
in  the  numerator,  3  (N-l)  in  the  denominator  where 
.is  the  number  of  elements  in  the  array.  The  figure 
is  for  an  array  of  J  3  e  1  ements-.  ^ 

False  alarms  per  day  as  a  function  of  probability 
of  detection  for  certain  fixed  signal  to  noise  ratio, 
on  the  array  beam,  and  for  3  degrees  of  freedom  in  the 

J™fratof’  J  (N_1)  ln  the  denominator  where  N  is  the 
number  of  elements  in  the  array.  The  figure  is  for 
an  array  of  19  elements.  L 
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Figure  No. 


I  also  a  1  aims  per  clay  as  a  1  unction  of  probability 
o!  detection  for  certain  fixed  signal  to  noise 
latios  on  the  array  beam,  and  for  3  degrees  of 
freedom  in  the  numerator  3  (N-l)  in  the  denominator 
where  X  is  the  number  of  elements  in  the  array.  The 
f  igure  is  f  or  an  array  of  31  elements. 

Signal  to  noise  ratio  required  at  the  seismometer 
ior  90b  probability  of  detection  as  a  function  of 
the  number  of  seismometers  in  the  array,  for 
different  values  of  false  alarms  per  day. 

bandpass  prefilter  for  the  Fisher  detector. 

TFO  37-element  array.  Elements  3,  5,  7,  8,  12,  and 
lb  not  used  in  the  Fisher  detector. 

Detector  response  on  0732  April  25,  1970,  Trace  9 
shows  the  Kuril  beam  through  the  filter  shown  in 
1 igure  9.  Irace  10  is  through  a  broader  filter, 
and  trice  11  is  the  detector  operating  on  trace  9. 

Detector  response  on  0111  April  9,  1970.  Trace  9 
shows  the  Kuril  beam  through  the  filter  shown  in 
Figure  9.  Trace  10  is  through  a  broader  filter, 
and  trace  11  is  the  detector  operating  on  trace  9. 

Detector  response  on  0644  April  25,  1970.  Trace  9 
shows  the  Kuril  beam  through  the  filter  shown 
in  I igure  9.  Trace  10  is  through  a  broader  filter, 
and  trace  11  is  the  detector  operating  on  trace  9.' 

Detector  response  on  0052  April  8,  1970,  showing 
breakdown  of  detector  for  large  events.  Trace  9 
shows  the  Kuril  beam  through  the  filter  shown 
m  Figure  9.  Trace  10  is  through  a  broader  filter, 
and  trace  11  is  the  detector  operating  on  trace  9.' 
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f  igure  Title 


Figure  So. 


Detector  response  on  1343  April  19,  1970  showing 
detector  response  to  off-beam  event.  Trace  9 
shows  the  Kuril  beam  through  the  filter  shown 
in  figure  9.  Trace  10  is  through  a  broader 

haJeT.  traCC  U  15  the  de"--ctor  operating  on 

3  db  Contour  of  short-period  TFO  Kuril  beam 
Solid  dots  show  o.9  magnitude  events  with  detector 

values. 8rCater  than  /*°*  °Pen  dots  have  lesser 

Piobabilitv  of  detection  as  a  function  of  magnitude 

of8?)!1  tn  n°1Se’  Peak  t0  Peak  ratios  as  a  function 
of  the  corresponding  value  on  the  Fisher  detector 
Arrow  indicates  F  =  10.3  threshold. 

Probability  that  detector  trace  exceeds  F  as  a 
function  of  T .  Theoretical  curves  for  A  =  0  0 
and  1.0;  experimental  curves  for  three  different 
days  as  noted,  for  each  day  histograms  are  given 
for  successive  two-hour  periods. 

Probability  that  F  >  2.0  on  the  Kuril  beam  as  a 
unction  of  time  for  the  period  August,  1970. 

Probability  that  F  >  2.0  as  a  function  of  the 
average  of  the  six  highest  velocity  flags  from 
United  States  weather  maps  inside  the  3  db  TFO 
contour.  Letters  indicate  points  for  which  maps 
are  presented  in  Figure  22.  Solid  points  are  from 
August  16  through  24.  Other  points  as  indicated. 
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figure  I"  i  tie 


dinted  States  Heather 
indicated  bv  letter  a 
August  16,  1 9 ” 0 . 


Bui eau  map  for  the  points 
in  Figures  19  and  20, 


United  States  Heather 
indicated  by  letter  b 
August  19,  1970. 


Bureau  map  for  the  points 
in  Figures  19  and  20, 


United  States  Heather 
indicated  bv  letter  c 
August  24,  1970. 


Bureau  map  for  the  points 
in  ligures  19  and  20, 


Figure  Ho. 
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,;!inoSn!rahty  parameters ,  for  a  31-element 
and  90  degrees  of  freedom.  element 
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INTRODUCTION 


The  purpose  of  this  report  is  to  describe  the  background 
and  the  present  status  of  the  Fisher  detector  of  seismic  signals 
and  to  compare  its  efficiency  with  other  types  of  detectors. 

Two  basic  types  of  automatic  detectors  of  seismic  signals 
have  been  investigated  in  the  past.  In  the  first,  a  short  term 
average  of  the  rectified  or  squared  seismogram  is  compared  to  a 
longer  term  average;  when  the  ratio  exceeds  a  threshold,  a  detec-  - 
ti°n  is  declared.  In  the  second  type,  a  checking  procedure  is 
either  explicitly  or  implicitly  programmed  to  see  whether  several 
independent  seismometers  give  a  detection  within  specified  time 
windows.  This  checking  procedure  is  analogous  to  the  procedures' 

of  a  human  operator  and  can  reject  data  spikes,  local  events,  and 
off-beam  events. 

The  detector  implemented  by  IBM  personnel  at  SAAC  (1969)  is 
an  example  of  the  first  type  of  detector.  A  full  LASA  beam  is 
rectified  and  the  result  is  averaged  over  approximately  2  seconds. 
This  is  output  to  a  long-term  exponential  smoothing  algorithm  in 
which  the  n  short-term  exponential  average  is  equal  to  (1-a) 
times  the  (n-l)st  value  plus  a  times  the  current  rectified  average. 
The  time  constant  is  typically  set  at  60  seconds.  A  fixed  value 
of  the  ratio  of  the  short-term  average  to  the  exponential  average 
is  the  threshold  for  declaring  a  detection. 

The  IBM  detector  can  be  triggered  by  off-beam  and  local  events. 
It  is  protected  from  triggering  on  data  spikes  because  the  raw  data 
is  limited  to  a  fixed  equivalent  ground  motion  (11  my)  and  when 
such  a  value  is  beamed  and  averaged  it  does  not  create  excursions 
large  enough  to  cause  a  detection. 

An  earlier  detector,  developed  by  Lincoln  Lab  personnel  (1965), 
calculated  a  similar  shorthand  long-term  average  on  the  traces 
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ot  the  center  seismometers  of  the  E  and  F  r^ng  subarrays  of  LASA. 

If  a  threshold  crossing  occurred  on  several  of  the  eight  channels 
within  20  seconds  of  one  another,  a  detection  was  declared.  This 
detector  will  adequately  reject  data  spikes,  but,  of  course  triggers 
on  al1  seismic  events  regardless  of  azimuth.  By  narrowing  the 
arrival  windows  to  suit  certain  move-out  velocities,  one  would 
trigger  detections  only  on  arrivals  from  a  fixed  epicentral  region. 
However,  this  refinement  was  not  programmed  by  Lincoln  Labs.  The 
other  disadvantage  of  this  particular  detector  is  that  it  makes  no 
use  of  the  improvement  in  signal  to  noise  on  the  beam. 


The  Fisher  detector  was  first  suggested  by  Melton  and  Bailey 
(1957).  The  essence  of  this  detector  is  that  the  power  on  the 
beam  is  divided  by  the  power  of  the  difference  between  the  beam 
and  the  individual  channel  seismograms  averaged  over  all  channels, 
lo  express  this  mathematically,  consider  an  infinite  velocity  signal, 
then  let  the  individual  traces  be  given  by  Si  (t)  .  In  general  the 
traces  will  be  shifted  to  form  a  beam  given  by  S  =  (1/N)  £  S-.  The 
residual  error  on  each  channel  is  given  by  S.  -  S  where  * 

has  been  appropriately  shifted  for  non-infinite  velocity  beams. 

The  expression  for  the  Fisher  detector  on  N  channels  is  given  by 


F 


2(N--)  *  £(SrS)2 


(N-l) 


(lsi) 2 

Z"(S.-S)2 


(1) 


where  the  bar  indicates  that  a  function  has  been  averaged  over  some 
time  window,  i.e. 


S. 
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I  ilc  express  i  on  (1.1) 


can  be  transformed  to 
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UhiCh  is  Usually  simPler  compute.  If  the  individual  channels 
have  non- zero  means,  a  term  must  be  added  to  the  numerator  and  the 
result  is 


N-  1 
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(Ssp 


f[(ZSj)2J 


2  s; 
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I'hon  a  signal  arrives,  the  output  of  tins  detector  increase, 
lor  two  reasons,  hirst,  the  bear,  power  in  the  numerator  increase, 
■econd,  lor  high ■ signal- to-noiso  ratios,  the  denominator  nearlv 
unices.  Note  that  this  Kill  not  happen  unless  the  signal  ha,  the 
proper  move-out ;  otherwise  the  shifted  beam  cannot  cancel  the 
individual  seismometer  trace  so  the  denominator  remains  large 
In  this  way  the  Fisher  detector  checks  that  the  signal  is  recorded 
b,  several  seismometers.  For  an  off-beam  signal  both  the  numerator 

and  denominator  increase  proportionally  and  there  is  no  increase 
m  the  ratio. 

As  discussed  in  the  next  section,  these  useful  qualitative 
properties  are  complemented  by  the  fact  that  the  ratio  ha,  a  well 
studied  statistical  distribution.  The  ratio  follows  the  central  F 
l st r i but  ion  in  the  absence  of  signal,  and  the  non-central  F  distri- 
bution  in  the  presence  of  signal.  Thus  it  is  possible  to  p.odict 

7tabllS"°d  tl,0°ry  th°  l,robabilit>'  of  detection  and  false  alarm 
rate  for  any  combination  of  bandwidth,  signal  detection  window  and  ' 
number  ol  seismometers. 


Although  lelton  and  Bailey  (1957)  first  discussed  the  possi¬ 
bility  of  using  the  I isher  criterion  as  a  detector,  they  did  not 
appl\  it  to  seismic  data.  Booker  (1965)  was  the  next  to  study  the 
fisher  detector.  His  work  was  implemented  in  a  program  on  a  CDC  1604. 
He  showed  that  the  Fisher  detector  had  considerable  promise,  and 
that  it  was  superior  to  a  polarity  coincidence  detector  proposed  by 
Melton  and  Karr  (1957).  However,  Booker  had  available  only  data 
from  closely-spaced  observatory  seismometers.  Noise  recorded  by 
these  seismometers  is  correlated  in  the  signal  band,  which  is 
contrary  to  the  hypotheses  of  the  Fisher  detector.  Booker  recognized 
this  fact,  and  saw  the  effects  in  his  results.  He  expressed  the  hope 
that  it  would  be  possible  to  allow  for  the  correlated  noise  by  finding 
and  setting  a  higher  threshold  than  would  be  predicted  from  theory, 
and  that  this  threshold  would  be  stable. 

An  auxiliary  processor  was  built  by  Texas  Instruments  to  imple¬ 
ment  the  Fisher'  detector  in  real  time,  and  an  evaluation  of  it  was 
performed  on  line  at  CPO  by  Edwards,  Benno,  and  Creasey  (1967).  They 
found  that  no  stable  threshold  could  be  set,  and  advocated  that  a 
modification  be  made  to  the  processor  to  allow  for  an  adaptive 
threshold.  From  our  present  vantage  point  we  can  see  that  correlated 
noise  between  the  nineteen  closely  spaced  CPO  seismometers  was  the 
cause  of  this  difficulty,  and  that  the  detector  would  have  been 
stable  if  only  the  three  most  widely  spaced  seismometers  with  a 
spacing  of  2.9  kilometers  had  been  used.  The  closest  spacing  for 
the  full  array  was  0.3  kilometers. 

shumway  and  Husted  (1969)  showed  how  to  extend  the  results  of 
Booker  to  colored  spectra,  and  also  showed  how  to  compute  the 
probability  of  detection  given  the  signal- to-noise  ratio. 

In  March  of  1969  the  TI  auxiliary  processor  was  installed  on 
the  37-element  TFO  array  which  has  a  minimum  spacing  of  4  km.  For 
the  first  time  it  was  operating  on  data  which  was  uncorrelated  in 
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the  signal  band.  The  principal  conclusion  of  this  report  is  that 
at  TFO  the  Fisher  detector  performs  in  agreement  with  theory,  and 
is  an  excellent  detector  when  installed  at  an  array  whose  seismo¬ 
meters  are  spaced  widely  enough. 

If  two  signals  arrive  at  the  same  time,  one  of  them  must  be 
regarded  as  noise  and  detection  is  more  difficult.  A  useful 
extension  of  the  Fisher  detector  would  therefore  be  detection  of 
two  or  more  signals  simultaneously  present.  Shumway  and  Dean  (1968) 
have  presented  the  theory  required  to  produce  minimum  variance, 
undis torted  estimates  for  each  of  a  multiple  set  of  signals  imbedded 
m  a  collection  of  independent  stationarily  correlated  noise  series. 
For  the  case  of  Gaussian  noise,  Shumway  (1970)  has  presented  the 
asymptotic  maximum  likelihood  detector  for  each  of  the  multiple 
signals  in  terms  of  the  classical  F  statistic.  Work  is  presently 
underway  to  develop  the  probability  of  detection  and  false  alarm 
statistics  for  two  signals  in  a  horizontal  array:  and  to  test  the 
theory  on  real  data.  This  detectof  could  be  used  to  detect  explo¬ 
sions  set  off  on  purpose  during  a  large  earthquake,  and  to  detect 
and  estimate  the  long  period  waves  from  two  events  for  discrimi¬ 
nation  purposes.  Long  period  waves  frequently  overlap  in  time  since 
they  are  of  long  duration.  The  estimators  are  not  the  same  as  the 
so-called  multichannel  filter  which  treats  one  of  the  signals  as 
noise  and  estimates  a  noise  correlation  or  spectral  matrix  to  a^ply 
in  the  single-signal,  weighted  minimum  variance,  undistorted 
procedure.  It  has  the  advantage  that  a  multichannel  filter  can 
only  be  applied  to  estimate  the  second  event  and  then  only  if  a 
substantial  portion  of  the  first  event,  substantially  uncontaminated  by 
the  second  event, is  available  to  serve  as  the  noise  sample.  The 
estimators  developed  by  Kobayashi  and  Welsh  (1969),  called  by  them 

optimum  fidelity  estimators,  are  special  cases  of  those  developed 
by  Shumway  and  Dean. 
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THEORY 


The  basic  derivations  of  the  Fisher  detector  can  be  found 
in  the  references  cited  in  the  Introduction.  Here  we  shall  briefly 
review  the  notation,  and  present  some  new  theoretical  results  on 
the  detection  capability  of  arrays  using  the  Fisher  detector.  This 
theory  is  perhaps  of  more  than  academic  interest  since  the  evalua¬ 
tion  of  the  Fisher  detector  in  the  next  section  shows  it  to  behave 
in  good  agreement  with  theory  and  to  be  an  excellent  detector  which 
might  well  be  implemented  in  real  systems. 

The  basic  hypotheses  under  which  the  Fisher  detector  theory  is 
valid  are: 

1.  Identical,  deterministic  signal  on  each  channel. 

2.  Independent,  identically  distributed,  stationarily  auto- 
correlated,  Gaussian  noise  series  on  each  channel. 

The  detector  is  actually  only  sensitive  to  the  assumptions  of 
equal  signal  waveform  and  of  uncorrelated  noise  between  channels.  Non- 
Gaussian  noise  is  permissible  to  the  degree  that  the  central  limit 
theorem  guarantees  that  the  output  of  a  narrow  band  filter  operating 
on  any  distribution  tends  to  a  Gaussian  distribution.  The  noise  need 
be  stationary  only  over  the  time  window  surrounding  the  signal, 
which,  in  the  present  application  is  only  three  seconds.  The  equality 
of  the  noise  autocorrelations  in  the  detection  band  is  needed  in 
order  that  power  in  the  residual  noise  will  be  the  same  on  each 
channel  so  that  the  average  residual  noise  power  will  follow  on 
chi -  squared  distribution.  If  this  assumption  is  violated  the  detector 
may  not  follow  the  predictions  of  theory,  but  a  stable  threshold 
can  probably  be  set. 

Under  the  assumption  that  no  signal  is  present,  the  probability 
of  the  detection  trace  exceeding  a  given  value  of  F,  say  F  ,  is 
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f 


expressed  in¬ 


i’ 


1 


(4) 


l.qu.ition  (4)  gives  the  probability  that  anv  value  of  [•  drawn 
rro»,  an  I  distribution  with  X,  and  X,  degrees  of  freedom,  and  with 
a  non- con t rai i ty  parameter  (a  function  of  the  signal  to  noise  ratio, 

°  tcr#>  oxcccJ  Fo-  fo*  i»  reference  to  equations  12)  and  (4,  : 

■s  the  number  of  degrees  of  freedom  in  the  numerator,  X,  is  the 
degrees  of  1  reedom  for  the  denominator,  and  1  equal  to  0  in  this 
case,  is  the  non-centrality  parameter. 

Ihe  symbol  for  the  non- central i ty  parameter  is  chosen  to  be  ,n 
agreement  with  that  in  Abramowitz  and  Stegun  (1964).  The  relation 
between  a  and  the  signal  to  noise  ratio  on  the  beam  is  given  by 

A  = .  'IJTfS/NJ  *  S humway  and  llusted,  and  Scheffe  (1959)  define 
their  non -centrality  parameter  as  62 
and  61"  is  A  =  6  . 


and  the  relation  between  A 


Tho  degrees  of  freedom  for  X,  for  the  numerator  may  be  calcu¬ 
late, 1  as  follows:  samples  per  second  x  (band  width/folding  frequenevj 
x  signal  time  window  -  2  x  band  width  x  signal  time  window  .  2BT. 

Ihe  degrees  of  freedom  Nj  for  the  denominator  is  given  by  fx-lj  times 
Nl>  where  \  is  the  number  of  channels. 

The  false  alarm  rate  in  false  alarms  per  day,  with  a  threshold 
set  at  lo  and  a  signal  time  window  of  length  T  seconds  is  given  by 


False  Alarms  per  day  =  P(F \’2,  0]  >  r  }  x  1 


4x60x60 


P  A 

-  x  8.65  x  104  . 


f5) 
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i  equation  may  unde  res t  imate  the  lalse  alarms  since  it  assumes 
tli‘t  ii  neither  oi  two  adjacent  windows  of  length  x  give  a  false 
ilar;  ,  then  neither  will  an  intermediate  window  composed  of  onc-half 
inch.  However,  considering  that  this  estimate  would  he  highly 
i.  i  .  t.  1  a  t e d  with  estimates  in  the  other  two  windows,  the  increase 
a  t  e  lalse  ulaim  rate  would  certainly  he  much  less  than  a  factor 
lKn-  1  '  1  il(-‘  lalse  alarm  rate  fl-ARj  corresponding  to  F  calculated 
I  i  oi,i  equation  (SJ  there  corresponds  a  probability  of  detection  for  a 
l  i.wd  s  i  gna  1  -  t  o  -  no  i  s  e  level.  This  probability  is  given  by 


'■{'■[•V  'a-  1  f-S/N)J  ■  !,.( 

where  IS/N)  is  given  hv  A  =  2BT(S/.\J2 

fo  calculate  actual  probability  of  detection- fa lse  alarm 
curves  it  is  obviously  necessary  to  calculate  the  central  and 
non-central  f  distributions.  Tables  of  these  distributions  exist, 
however;  the  smallest  value  of  F  generally  tabulated  is  0.01,  and 
1 rom  (5)  with  a  time  window  of  5.0  seconds,  this  would  imply  255 
lal  e  alarms  per  day.  Thus  we  need  to  evaluate  these  distributions 
lor  ranges  outside  the  customary  tables.  The  methods  for  doing  so 
are  given  in  Appendix  1. 

Figure  1  shows  the  P  value  corresponding  to  a  particular  false 
1 1  a t a  rate  ior  a  succession  oi  numbers  of  seismometers  in  an  array. 

I  l  example,  with  51  elements  in  the  array,  and  a  false  alarm  rate 
of  0.1  per  day,  the  threshold  value  for  1:  i  s  10. 92,  using  the  left- 
hand  scale.  This  figure,  and  all  following  ones  assume  2 BT= 5 . 0  degrees 
of  freedom  in  the  numerator  (see  first  paragraph  in  Section  5J  and 
5(.\-l)  in  the  denominator  i.e.  Nj  =  5,  N  =  ( N  - 1 J  5 . 

I  igures  2  through  7  show  the  false  alarm  rate  as  a  function  of 
the  probability  of  detection,  bach  figure  is  for  a  different  size 
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array;  and,  on  each  figure,  each  curve  is  for  a  different  signal- 
to-nojsc-  ratio  on  the  beam,  for  example,  figure  7  is  for  31  e  lcments 
and,  with  a  s  igna  1  -  to- no i se  ratio  of  3.0  on  the  beam,  one  has  a 
90  percent  probability  of  detection  with  a  false  alarm  rate  of  0.22 
per  day.  he  may  mention  here  that  the  si  nal-to-noisc  ratio  is 
defined  as  the  ratio  of  the  square-root  of  the  signal  power  divided 
I’)'  the  square  root  d  the  noise  power  in  the  frequency  band  and 
time  window  under  consideration.  If  we  assume  that  the  rins  noise 
in  the  signal  window  is  equal  to  the  rms  noise  before  the  signal, 
and  that  the  signal  itself  is  typically  a  segment  of  a  sine  wave; 
then  to  translate  the  s ignal - to- no i s e  values  given  in  this  report 
to  values  in  terms  of  zero-peak  signal  divided  by  rms  noise  one 
would  multiply  by  1.414.  To  the  extent  that  the  signal  cannot  be 
regarded  as  a  sine  wave  in  the  3-second  window  used  in  this  report 
the  multiplying  factor  would  be  greater  to  allow  for  the  greater 
amplitude  of  one  out  of  the  three  cycles.  It  is  generally  agreed 
(e.g,  Clark,  1968)  that  the  conversion  between  peak- to  -  peak/ 2  rms 
s i g na 1 - to  -  no  i  se  values  and  peak-to-peak  signal/maximum  peak-to-pcak 
noise  is  performed  by  dividing  the  former  values  by  3.0.  Thus,  to 
convert  s ignal- to-noise  ratios  used  in  this  report  to  those  measured 
by  the  ratio  of  peak-to-peak  values  of  signal  and  noise  one  would 
multiply  by  1.414/3.0  equal  to  0.47  ~  0.5. 

figure  8  shows  the  s igna 1  -  to- noi se  ratio  required  at  the 
seismometer  level  for  90  percent  probability  of  detection  as  a 
function  of  the  number  of  instruments  in  the  array.  There  arc  four 
curves  for  four  different  false  alarm  rates,  he  see  for  example  that 
for  51  instruments  and  for  a  false  alarm  rate  of  0.1  per  day,  that 
the  s igna 1- to-noise  ratio  at  the  seismometer  level  must  be  0.51.  Also 
in  Figure  8  is  an  indication  of  the  range  in  magnitude  and  dB 
corresponding  to  a  range  in  s igna 1  - to-noise  ratios.  he  note  that 
for  all  arrays  with  more  than  20  elements  that  an  increase  in  false 


iliin  rate  of  a  factor  of  1,000  yields  a  decreases  in  the  magnitude 
detected  of'  only  0.13. 


The  s ignal - to-no ise  required  for  a  fixed  probability  of  detec¬ 
tion  as  a  function  of  the  number  of  seismometers  is  proportional 
to  X  -  for  N  greater  than  20.  Although  there  is  probably  a  simple 
explanation  for  this  fact,  we  have  not  found  it  as  yet.  It  is  clear 
f  rom  the  t igure  that  the  s ignal- to-noise  ratio  required  decreases 
much  more  rapidly  than  N  for  small  values  of  N.  This  reflects 

the  decreasing  degrees  of  freedom  in  the  denominator  of  the  1 
statistic.  A  power  detector,  for  which  the  noise  estimate  is 
determined  by  a  long  window  on  the  beam  in  front  of  the  signal, 
would  not  depart  from  N1/2  in  this  way  and  therefore  may  have 
advantages  for  small  arrays. 


/ 
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PROCEDURES 


fi8“"  9  sho"s  the  'tCr-Fishcr  pref i  1  ter  used  at  TFO.  It  is  a 
I>oint  symmetric  filter  designed  using  program  P  [  L  PLOT  2  with 
the  parameters  as  indicated  on  the  figure.  It  is  interesting  that 
to  force  a  null  at  0.2  Hz  to  cancel  the  microsetsms,  the  parameter 
wluch  was  adjusted  was  the  number  of  points  in  the  filter  lor  a 
time  Window  of  3  seconds,  the  number  of  degrees  of  freedom  allowed 
7  this  filter  is  3.6,  assuming  a  white  signal  spectrum  in  the 
pass-band.  In  0ur  calculations  we  have  rounded  3.6  to  3.0  because 
t  e  algorithms  for  calculating  the  F  distribution  do  not  accept 
non-integral  values  of  iq  or  iq  and  because  the  signal  and  noise 
spectra  are  not  truly  white. 

It  is  important  to  determine  properly  the  length  of  the  signal 
time  window.  Intuition  suggests  that  the  window  should  be  equal  in 
length  to  the  visible  portion  of  the  P-wave.  Booker  (1965 j  reported 
that  he  found  this  to  be  true.  Edwards  et  al  (1967)  reported  that 
the  typical  small  signal  was  3  seconds  in  duration;  and  that  the 
detector  seemed  unstable  when  a  window  of  only  2  seconds  was  used. 

■  ri  the  present  study  we  used  a  time  window  of  3  Seconds.  Booker 
n5°  presented  results  to  show  that  a  samplxng  rate  of  10  sps  was 
adequate ,  however  at  TFO  the  processor  operates  at  20  sps.  Thus  the 
time  window  for  the  auxiliary  processor  was  set  to  be  60  points. 

A  full  presentation  of  the  MCF  and  auxiliary  processor  control 
data  for  one  set  of  MCF  beams,  including  a  Fisher  detector  on  a 
Kuril  beam  is  given  in  Appendix  II,  together  with  a  history  of 
the  Fisher  beams  recorded  on  film  and  tape  up  to  the  present  time, 
so  that  the  interested  researcher  can  use  the  data  easily. 

figure  10  shows  a  map  of  the  37  element  TFO  array  used.  Elements 
3,  5,  7,  8,  12,  and  16  are  not  in  the  beam  used  by  the  Fisher 
detector,  leaving  a  total  array  of  only  31  elements. 
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The  output  ot  the  Fisher  detector  is  recorded  both  on  magnetic 
tape  in  digital  form,  and  on  develocorder  film.  In  measuring  ampli¬ 
tudes  from  the  develocorder  it  is  important  to  follow  systematic 
procedures  because  even  an  error  of  0.1  cm,  as  seen  on  the  film 
viewer,  is  significant.  Since  the  film  can  slip  up  and  down  in 
the  viewer,  and  since  different  viewers  can  yield  slightly  different 
measurements  for  the  same  point,  it  is  essential  always  to  use  the 
same  viewer,  to  allow  the  film  to  move  only  in  one  direction  through 
the  viewer,  and  to  make  all  measurements  along  the  vertical  center 
line  of  the  viewer. 

Calibrations  on  the  film  were  maintained  by  finding  the  lowest 
value  of  the  Fisher  in  a  30  minute  period,  three  times  a  day;  2-2:30, 
12-12:30,  and  22-22:30.  Examination  of  the  digital  tape  shows  that 
the  Fisher  goes  to  zero  at  least  every  10  minutes,  and  so  these 
minimum  values  should  be  at  the  same  place  on  the  viewing  screen  if 
the  calibration  is  to  be  held  from  beginning  to  end  of  the  film. 
Readings  are  not  accepted  unless  the  minima  are  within  0.1  cm  of  one 
another.  Conversion  of  the  Fisher  readings  (in  centimeters)  are  con¬ 
verted  to  F  values  by  creating  a  full  scale  output  from  the  Fisher 
detector  which  is  registered  on  the  film,  and  on  the  magnetic  tape. 
Calibrations  performed  on  30  June  and  3  August  1970  were  within 
0.1  cm  of  one  another  on  the  film  and  within  four  counts  out  of  715 
on  the  digital  tape.  The  Fisher  output  which  is  digital,  is  con¬ 
verted  to  analog  form,  transmitted  to  the  Astrodata  equipment, 
redigitized  and  recorded. 

An  unfortunate  feature  of  the  Fisher  detector  as  implemented  is 
that  overflows  occur  for  large  events,  either  on  or  off  beam.  This 
is  manifested  on  the  develocorder  film  by  a  discontinuity  in  the 
Fisher  trace.  The  same  feature  was  noted  by  Edwards  et  al  (1967); 
they  state:  "Th^  MCF  Auxiliary  Processor  system  is  limited  to  a 
12-bit  (66  db)  dynamic  range  on  input...  Large  signals  of  interest 
were  clipped  on  input  or  during  intermediate  computations..." 
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Figure  11  shows  the  detector  response  at  0732  April -25,  1970. 
Traces  0-3  are  beamed  for  a  A  of  30°,  and  azimuths  of  0,  90,  180, 
and  270  degrees.  Traces  4-7  are  for  a  A  of  100°  and  the  same  four 
azimuths.  Trace  8  is  for  infinite  velocity  and  trace  9  is  aimed 
at  the  Kuril  Islands,  a  A  of  70°  and  an  azimuth  of  315°.  All  of 
these  beams  are  filtered  with  the  same  23  point  filter  shown  in 
Figure  9.  Trace  10  is  the  Kuril  beam  through  a  broader  filter. 

Trace  11  is  the  Fisher  detector  itself,  and  trace  .12  is  the 
infinite  velocity  beam  of  the  crossed  array  at  TFO.  We  note  that 
the  signal  is  clearly  detectable  on  traces  9  and  10,  and  that  the 
detection  trace  rises  smoothly  to  a  large  value.  Figure  12  shows 
the  arrival  of  another  we  1 1 - recorded  event  at  0111  on  9  April  1970. 
Data  soikes  are  responsible  for  the  patterns  of  "arrival.  "  to  the 
left  and  right  of  the  true  arrival.  Note  that  the  detector  is 
unaffected.  Figure  13  shows  the  behavior  of  the  detector  on  a  weak 
event;  and  Figure  14  shows  how  it  breaks  down  due  to  hardware  over¬ 
flows  for  a  large  event  on  beam.  Figure  15  shows  the  behavior  of  the 
detector  during  a  large  arrival  on  another  beam.  We  note  that  the 
trace  flattens  out,  suggesting,  as  is  true,  that  it  would  be  more 
difficult  to  detect  an  on-beam  signal  during  this  time  interval.  We 
note  also  the  small  sawtooth- 1 ike  pattern  on  the  Fisher  trace  indi¬ 
cating  over-flow  conditions  during  the  largest  arrivals. 


results 


r,e„rc  16  shows  the  3  dB  contour  of  the  TFO  short  period 
or.l  heam.  Events  used  were  those  reported  during  the  periods 

,',57™  ,  1-15  J“"e  197°’  fr0"  thc  ro«i<>n  *>°™ded  V 

-0 WON  and  110E-160W.  The  arrival  time  of  each  event  at  the 

center  element  was  predicted,  and  the  maximum  value  of  the 

IS  er  in  a  20-second  window  centered  at  this  time  was  recorded. 

Assuming  a  threshold  nt  E  -  in  i  j 

t  ,,  "  at  F  -  10-3  and  assuming  that  all  breakdown 

points  would  have  gone  over  this  threshold  if  there  had  been 
no  overflow  one  may  generate  Figure  17  which  gives  the  prob- 

y  °  ctectl°n  as  a  function  of  magnitude.  This  figure 
suggests  that  the  incremental  90  percent  detection  threshold 
magnitude  is  4.15  ±  .05. 

An  •  value  of  10.3,  by  Figure  1  for  n  - 

’  y  tor  N  -  31,  corresponds  to 

coretical  false  alarm  rate  of  0.2  per  day.  The  develocorder 
s  or  13  days  selected  at  random  between  08  April  and  15 
une  were  scanned  for  all  excursions  of  the  Fisher  trace  greater 

.  =  10-3'  F°r  10  °f  th0  ^  ‘hare  were  no  false  alarms,  and 

here  was  one  each  for  three  other  days  for  a  false  alarm  rate 

0.2.1  per  day.  The  expected  number  of  detections  that  are  true 
events  not  reported  can  be  calculated  as  approximately  .06  per 
day.  There  were  seven  excursions  over  F  •  10.3  during  the  13  days 
-which  were  traced  to  true  events  in  either  the  Philippines  or 
southern  Alaska;  outside  of  the  3  dB  contour. 
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Thus,  in  summary,  there  was  an  average  of  2.0  true  detections 
per  day  with  a  90  percent  threshold  magnitude  of  4.15  +  .05,  0  23 
apparent  false  alarms  per  day,  and  0.54  off-beam  detections  per 
day.  Of  the  apparent  false  alarms,  .06  might  he  unreported  events, 
leaving  an  estimated  true  false  alarm  rate  of  0.17. 

figure  18  shows  the  s  ignal - to-no ise  value  on  the  Kuril  beam  as 
a  function  of  the  Fisher  value  for  all  events  with  Fisher  values 
greater  than  10.3,  and  a  few  day's  data  for  smaller  values.  The 
s i gna 1 - lo -noise  value  was  determined  from  film  measurements  as  the 
ratio  of  the  peak-to-peak  signal  amplitude  to  the  peak-to-peak  noise 
amplitude  in  the  previous  60  seconds.  The  discussion  in  the  Theory 
section  shows  that  we  should  multiply  these  signal- to-noise  values 
by  a  factor  close  to  2.0  to  convert  them  to  the  signal-to-noi so  values 
used  elsewhere  in  this  study.  This  would  imply  that  the  10.5  thres¬ 
hold  corresponds  approximately  to  a  signal-to-noise  rms  power  ratio 
of  4.0  on  the  beam.  Fxamination  of  Figure  7  shows  that,  for  31 
elements  a  false  alarm  rate  on  the  order  of  0.2,  corresponding  to 
a  Fisher  threshold  of  10.3,  with  90  percent  probability  of  detection, 
does  indeed  require  a  rms  signal-to-noise  value  of  approximately  3.0. 
Thus  we  see  that  the  detector  behaves  as  predicted  from  theory  when 
examined  from  this  alternate  direction. 

Figure  19  is  a  plot  of  the  probability  that  the  Fisher  trace 
exceeds  a  certain  F  value,  as  a  function  of  F.  The  two  solid  curves 
are  theoretical  for  non-centrality  parameters  of  A  =  0  and  1.0. 

Also  on  the  curve  are  experimentally  determined  points,  each  from 
two  hours  of  data  computed  using  program  FISHDET2  from  F  values 
recorded  on  digital  tape  at  TFO .  We  see  that  the  histograms  from 
successive  two-hour  periods  of  data  are  quite  close  together,  showing 
that  the  detector  statistics  are  stable.  However,  over  longer  periods 
of  time,  the  statistics  can  be  quite  different.  These  differences 
are  significant  as  can  be  seen  by  the  vertical  error  bars  giving 
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the  99  percent  Ko  linogorov- Smi  rnov  confidence  interval.  If  any 
two  of  these  histograms  are  drawn  from  the  same  population,  then 
th>  distance  between  them  at  every  value  of  1'  will  be  less  than 
tli  is  interval  99  percent  of  the  time. 


he  observed  that  large  probabilities  correlate  with  storms 
in  the  Kurils  region.  Texas  Instruments  (1967)  have  also  found 
storms  to  be  a  source  of  P-wavc  noise.  Figure  20  presents  the 
piobability  that  F  >  2.0  for  the  period  16-24  August  1970.  United 
-tates  heather  bureau  maps  revealed  that  a  low  pressure  system 
moved  across  the  beam  during  this  interval,  and  that  its  winds 
U‘ere  "t  a  maximum  on  the  19th.  Figure  21  shows  the  correlation 
between  the  average  of  the  six  highest  velocity  readings  inside 
the  .>  db  contour  as  reported  on  the  heather  Bureau  maps  for  the 
^amc  time  period,  and  the  probability  that  F  >  2.0.  The  solid 
dots  are  for  the  period  16-24  August,  and  the’other  points  are 
as  noted.  Uthough  there  is  considerable  scatter,  the  trend  is 
detimto.  Figures  22a,  b,  and  c  show  the  Weather  Bureau  maps 
lor  the  points  marked  by  a  letter  in  Figure  19. 


The  theoretical  value  for  P(F)  at  F  =  2.0,  \  =  0  is  p  =  0.1. 

1  ^e  Ulct  that  for  somc  da>'s  <  0.1  cannot  be  explained  by  a  non- 
centrality  parameter  different  from  0.0.  It  is  our  present 
beliel  that,  values  less  than  0.1  arise  on  days  when  significant 
storm  energy  arrives  at  TFO  from  locations  other  than  Kurils.  We 
leel  that  this  biases  the  data  so  that  it  is  more  difficult  for 
an  accidental  fluctuation  to  look  like  a  signal  from  the  Kurils, 
leading  to  a  low  probability  for  high  Fisher  values.  An  illustration 
o I  tllls  for  a  Inrgc  off-beam  signal  was  seen  in  Figure  15. 

The  low  points  in  Figure  20  correspond  to  days  when  the  low-pressure 
system  had  moved  out  of  the  beam. 


Leaving  aside  these  times  when  there  is  an  unknown  lower  false 
alarm  rate  (and  some  correspondingly  lowe>'  probability  of  detection), 
the  range  in  false  alarms  per  day  implied  by  non-centrality  parameters 
ranging  from  0.0  to  that  corresponding  to  the  largest  observed 
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p robab i 1 i ty  value  displayed  in  Figure  20,  approximately  A  =  1.0, 
is  0.2  to  1.90  false  alarms  per  day.  If  one  adjusted  the  threshold 
value  of  F  to  maintain  a  constant  false  alarm  rate  of  0.2  per  day, 
the  range  of  F  values  would  lie  between  10.3  and  12.7;  and  the 
corresponding  range  in  s ignal - to- no i s e  at  the  seismometer  level 
would  be  0.51  to  0.57;  a  range  of  0.05  magnitude  units.  Thus,  by 
adaptively  adjusting  the  threshold  one  might  hope  to  detect  events 
0.05  smaller  on  days  for  which  there  were  no  storms  than  on  davs 
when  storms  were  present.  The  problem  of  how  to  predict  the  false 
alarm  rate  and  probability  of  detection  in  the  presence  of  weak 
signals  from  directions  other  than  the  target  area  is  presently 
under  investigation  at  SDL  by  means  of  simulation. 

Figure  23  shows  the  threshold  F  value  as  a  function  of  the 

resulting  false  alarm  rate  for  a  series  of  non- central i ty  para¬ 
meters  . 

On  Figure  16  some  3.9  magnitude  events  are  plotted.  Those  with 
F  values  greater  than  7.0  are  plotted  as  solid  circles.  There  is 
some  tendency  for  greater  detection  closer  to  TFO ,  in  agreement 
with  the  fact  that  the  amplitude  falls  off  slightxy,  somewhat  less 
than  0.1  magnitude  units  on  average,  with  distance  in  the  range 
45-85°.  however ,  the  result  is  not  clearcut,  as  would  be  expected 
from  known  variances  in  the  determination  of  magnitude. 

In  a  study  by  Clark  (1970),  using  the  same  31  sensor  TFO  Kuril 
beam  with  a  0.7  to  2.5  Hz  filter,  an  experienced  analyst  obtained 
90  Pcrccnt  detection  threshold  of  4.35  +  .1  mb .  Adequate  records 
"ere  not  kcPt  to  enable  an  estimate  of  the  false  alarm  rate,  however 
the  analyst  (Paul  DcFries,  has  estimated  a  rate  comparable  to  that 
of  the  Fisher  detector.  Thus  the  Fisher  detector  is  apparently 
comparable  to  or  slightly  better  than  an  experienced  analyst. 
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There  is  no  available  data  on  which  to  base  a  comparison  with 
a  power  detector.  However  it  is  clear  that  the  false  ala™  rate 

hale  bC0"  "UCh  Meher  at  ™  *»  to  local  and  other  off-beam 
cats  Ihere  were  also  many  events  each  day  where  relaxation  pulses 
T  om  the  seismometers  would  have  triggered  a  threshold  detector- 

a  limiting  of  the  seismometer  traces,  as  done  in  the  IBm' 

-system,  might  have  solved  this  problem. 


\ 
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DISCUSSION 


The  Fisher  detector  has  been  shown  to  give  a  90  percent 
incremental  threshold  magnitude  of  4.15  +  .05  when  operating 
on  a  TFO  31-element  Japan-Kurils  beam.  The  corresponding  true 
false  alarm  rate  was  found  to  be  0.17  per  day.  Questions  which 
now  arise  are:  will  the  detector  work  in  agreement  with  theory 
when  fewer  elements,  e.g.  7,  13,  or  19,  are  used?  How  well  will 
it  work  when  the  fewer  elements  are  closely  clustered  yielding 
a  much  broader  beam  and  thus  more  detections  on-beam  per  day? 

What  is  the  detection  threshold  for  closer  teleseismic  distances? 
Can  the  detection  threshold  be  lowered  by  widening  the  band-pass, 
especially  toward  the  high-frequency  end? 

A  study  should  also  be  made  of  the  computing  capacity  required, 
with  the  latest  available  equipment,  to  implement  the  detector  on 
a  suite  of  world-covering  beams  in  a  realistic  system  for  arrays 
of  7,  13,  and  19  elements. 

The  Fisher  detector  should  be  monitored  through  the  winter 
and  hurricane  seasons  to  see  if  there  are  a  few  days  per  year  in 

which  the  detection  threshold  is  lowered  a  more  substantial  amount 
by  storms. 

Work  is  now  in  progress  at  SDL  to  extend  the  detector  so  that 
it  will  work  in  the  presence  of  another  signal.  It  is  our  belief 
that  this  will  prove  to  be  the  first  "high-resolution"  technique 
capable  of  "splitting"  two  simultaneous  events  to  a  degree  superior 
to  that  of  conventional  beamforming. 

Work  is  also  in  progress  to  determine  the  beam  pattern  of  an 
array  on  which  the  Fisher  detector  is  operating.  This  beam  pattern 
will  be  expressed  in  terms  of  probability  of  detection;  and  of 
course  the  beam  pattern  will  be  different  for  different  signal-to- 
noise  ratios.  Despite  this  complication  the  probability  beam  pattern 
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IS  the  one  truly  relevant  to  the  surveillance  problem.  It  enables 
one  to  cover  the  earth  with  beams  guaranteeing  a  particular 
probability  of  detection  for  some  fixed  magnitude. 

Using  the  probability  of  detection  false  alarm  curves 
deduced  m  this  study,  it  will  also  be  possible  to  deduce  the 
corresponding  curves  for  a  network  of  stations. 
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APPENDIX  I 


CALCULATION  OF  THE  F- DISTRI BUT  ION 


Tables  of  the  F  distribution  are  available  in  Scheffe,  (1959). 
However ,  they  are  tabulated  only  for^values  of  P  between  f) ’.99  and 
0.01.  In  addition  they  are  printed  in  a  form  somewhat  difficult 
to  use.  Equations  for  the  F  distribution  are  given  by  Abramowitz 
and  Stegun,  (1964).  We  used  formulas  26.6.7  and  26.6.18.  When  we 
first  attempted  to  use  the  formulas,  the  results  did  not  agree 
with  the  tables  in  Scheffg.  Ronald  Finkler  of  Institute  of  Defense 
Analyses  (personal  communication)  pointed  out  that  for  the  equations 
to  be  correct,  what  appears  as  F'  on  the  right  hand  side  of  equation 
26.6.18  must  be  replaced'by  F'vj/fv  +  2  j  ) 

Then  the  correct  equatibns  are: 
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In  the  computer  calculation  use  is  madl  of  the  fact  that 
successive  terms .in  26.6.7  may  be  calculated  recursively.  This 
formula  is  in  the  inner  J,r.op  between  statements  40  and  70 
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APPENDIX  II 


OPERATING  SCHEDULES  AND  PARAMETERS 
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false  alarms  per  day  as  a  func on  of  probability  of 
etection  for  certain  fixed  signal  to  noise  ratios  on 
the  array  beam,  and  for  3  degrees  of  freedom  in  the 
numerator,  3  (X-J)  in  the  denominator  where  N  is  the 
number  of  elements  in  the  array.  The  figure  is  for  an 
array  -of  7  elements . 
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figure  8.  Signal  to  noise  ratio  required  at  the  seismometer  for 
90o  probability  of  detection  as  a  function  of  the 
number  of  seismometers  in  the  array,  for  different 
values  of  false  alarms  per  dav . 
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I;  i  gure  9.  Bandpass  [ire!  iltcr  tor  the  Pisher  detector. 
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Figure  11.  Detector  response  on  0732  April  25,  1970,  ,rare  9 
shows  the  Kuril  beam  through  the  filter  shown  in 
figure  9,  Trace  10  is  through  a  broader  filter,  and 
t'-ce  11  is  the  detector  operating  on  trace  9. 
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Figure  12.  Detector  response  on  0.11  April  9,  1970.  Trace  9 
shows  the  Kuril  beam  through  the  filter  shown  in 
Figure  9.  Trace  10  is  through  a  broader  filter,  and 
trace  11  is  the  detector  operating  on  trace  9. 
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Figure  18.  Signal  to  noise,  peak  to  peak  ratio 
the  corresponding  value  on  the  i ishi 
indicates  V  =  10.3  threshold. 
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